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I.  INTRODUCTION 


A  problem  of  long  standing  in  the  design  of  radar  systems  has  been 
that  of  detecting  the  presence  and  estimating  the  position  of  a  ship  at 
sea.  Contemporary  radar  systems  can  exploit  a  variety  of  techniques, 
such  as  frequency  agility,  doppler  filtering,  time  domain  processing, 
polarization  diversity,  etc.  to  enhance  the  effective  target  contribution 
to  the  total  return  signal.  However,  the  rational  design  of  such  systems 
must  depend  on  a  detailed  understanding  of  the  scattering  characteristics 
of  the  target  (i.e.,  the  ship),  the  clutter  background  (i.e.,  the  sea 
surface),  and  their  interactions.  Furthermore,  the  in-si tu  testing  of 
such  systems  is  expensive,  and  the  results  are  often  hard  to  interpret 
because  one  cannot  reconstruct  the  precise  encounter  geometry  (including 
the  height  distribution  of  the  sea  surface)  at  the  moment  of  measurement. 

Thus  for  a  number  of  reasons,  it  is  desirable  to  utilize  mathe¬ 
matical  modeling  of  the  extended  target  complex  represented  by  the  ship 
in  the  presence  of  the  sea  surface.  Such  a  model,  once  it  has  been 
properly  validated,  will  also  permit  the  economical  design  of  a  direct 
physical  model  of  the  encounter  geometry.  During  the  past  two  years, 
under  Contracts  N60530-75-C-0216  and  N60530-76-C-0186,  the  ElectroScience 
Laboratory  has; 

(i)  designed  a  one  dimensional  physical  model  of  the 
ocean  surface  [1], 

(ii)  developed  techniques  for  calculating  the  scattering 
from  such  a  surface  and  relating  it  to  the  scattering 
from  the  ocean  [1,2,3], 

(iii)  developed  simplified  ship-sea  models  in  which  the 
interaction  term  (the  major  contribution  to  the  total 
signal)  can  be  identified  and  calculated. 


It  is  the  purpose  of  this  report  to  describe  these  simplified  ship 
sea  models,  including  the  sea-surface  models,  in  detail;  to  review  the 
techniques  for  calculating  the  several  terms  in  the  total  radar  signature 
from  such  models;  and  to  compare  a  number  of  calculated  signatures  with 
measurements  carried  out  at  the  Encounter  Simulation  Laboratory  of  the 
Naval  Weapons  Center. 


II.  REVIEW  OF  PREVIOUS  WORK,  AND  GENERAL  APPROACH 
TO  SHIP-SEA  MODELS 

A  number  of  significant  steps  have  already  been  taken  towards  the 
development  of  a  complete  near-field  ship-sea  model.  The  first  phase 
was  to  design  a  physical  surface  which  had  the  same  scattering  mechanisms 
as  the  sea  surface,  in  order  to  provide  both  an  inexpensive  physical 


model  of  the  sea  surface,  and  to  validate  mathematical  models  for 
determining  the  scattered  fields.  The  bistatic  scattering  behavior 
of  the  sea  surface  alone  is  quite  complex,  but  is  now  understood 
fairly  well  in  terms  of  the  "composite"  model.  In  this  model,  the 
scattering  near  the  specular  direction  (i.e.,  near  normal  incidence  for 
back  scatter)  is  controlled  by  the  slope  distribution  of  the  large 
scale  structure  of  the  surface.  This  part  of  the  scattering  is  then 
calculated  by  the  physical  optics  model,  and  provides  an  explanation 
not  only  for  the  near  normal  back  scatter  but  also  for  the  dependence 
of  brightness  temperature  on  polarization  and  wind  speed  (at  least  at 
the  lower  wind  speeds).  For  off-normal  incidence  the  back-scattering 
mechanism  is  the  "Bragg-scatter"  from  capillary  and  short  wavelength 
components  of  the  ocean  surface.  This  part  of  the  scattering  is  cal¬ 
culated  by  perturbation  theory  and  when  combined  with  the  known  height 
spectrum  of  the  ocean  it  explains  the  angular  and  polarization  dependence 
of  the  radar  cross-sect  ion,  and  the  weak  dependence  on  electromagnetic 
wavelength  and  wind  speed.  Thus  the  problem  in  specifying  a  suitable 
physical  surface  model  is  to  select  one  which  will  have  a  height  spectrum 
which  satisfies  the  Bragg  condition  for  the  given  experimental  geometry 
(for  back  scatter  away  from  normal  incidence),  has  an  effective  mean 
square  slope  equivalent  to  the  desired  ocean  surface  (for  forward 
scatter  or  back  scatter  near  normal  incidence)  and.  would  he  compatible 
with  the  geometrical  and  other  constraints  imposed  by  the  practical 
limitations  of  the  Encounter  Simulation  Laboratory  of  the  Naval  Weapons 
Center. 


Figure  2-1.  Composite  sea  surface-ship  model.  T 
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The  final  design  of  the  surface  is  reported  in  Reference  [1],  and 
is  illustrated  in  figure  .’-1.  The  relation  between  the  structural 
form  of  the  surface  and  the  scattering  mechanisms  of  the  ocean  are 
discussed  in  more  detail  in  Reference  [1],  In  general,  however,  the 
large  undulations  represent  the  large  scale  (gravity  wave)  structure 
of  the  ocean  and  the  small  corrugations  provide  Bragg  scatter  from 
the  small  scale  structure.  Three  mathematical  models  (physical  optics, 
perturbation  theory,  moment  method)  were  developed  to  calculate  the 
scattering  from  this  composite  surface,  and  these  are  described  in 
detail  in  References  [1-B].  It  was  found,  as  expected,  that  the 
physical  optics  model  provided  an  accurate  prediction  of  the  scattered 
field  near  the  specular  direction.  Away  from  t tie  specular  direction, 
it  was  expected  that  the  scattered  fields  would  be  calculated  by  per¬ 
turbation  theory,  but  it  was  found  that  while  the  general  level,  and 
the  polarization  dependence  were  correct,  the  detailed  field  structure 
was  not  well  predicted.  In  order  to  achieve  better  agreement  with 
measured  data,  the  modified  moment  model  [ .'  ]  was  developed.  In  this 
model,  the  moment  method  is  used  to  estimate  the  exact  surface  currents 
over  a  single  period  or  cell  of  the  small  scale  structure,  and  the 
total  signature  computed  from  the  sum  of  these  individual  cell  cur¬ 
rents.  This  modified  moment  method  produced  signatures  in  excellent 
agreement  with  the  measured  back-scattering  signatures.  Furthermore,  a 
comparison  of  the  measured  and  calculated  radar  cross-section  of  the 
composite  surface  with  the  known  cross-section  of  the  ocean  showed  that 
both  of  the  known  ocean  scatter  mechanisms  were  appropriately  modelled. 

Thus  the  next  step  in  developing  a  complete  model  for  the  interaction  of 
radar  system  with  a  ship  at  sea  is  to  produce  a  model  for  the  ship 
which  will  be  complex  enough  to  exhibit  all  the  major  scattering  mech¬ 
anisms  of  the  actual  situation,  but  still  simple  enough  so  that  the 
scattered  fields  may  be  readily  computed. 

In  undertaking  the  development  of  these  simplified  ship-sea  models, 
two  separate  approaches  have  been  taken.  The  first  approach  considered 
a  model  which  could  easily  be  implemented  using  the  composite  surface 
already  constructed  at  the  Encounter  Simulation  Laboratory  (Figure  B-l). 
This  would  permit  actual  measurements  (subsequently  carried  out  by  Naval 
Weapons  Center)  of  typical  signatures,  which  could  be  compared  with 
calculated  signatures  and  thus  validate  the  theoretical  basis  of  the 
model.  This  model,  the  techniques  for  calculating  the  signatures,  and 
a  comparison  with  the  measured  signatures,  are  described  in  more  detail 
in  Sect  ion  IV . 

The  second  approach  wa«  the  development  of  a  model  more  repre¬ 
sentative  of  an  actual  encounter  geometry,  utilizing  a  two  dimensional 
sea  surface  (discussed  in  section  III),  and  a  simple,  but  three  dimensional 
ship  model.  This  ship-sea  model  is  discussed  in  section  V. 


III.  MODELS  FOR  THE  SEA  SURFACE 

Before  undertaking  the  construction  of  the  ship-sea  models  in 
their  entirety,  i t  wi 1 1  be  useful  to  develop  an  appropriate  model  for 
the  sea  surface.  Three  different  sea  surface  models  have  been  con¬ 
sidered,  and  are  described  briefly  below. 

( i )  Distribution  of  Specular  Points 

It  is  known  [4]  that  the  specular  component  of  scatter  from  the 
sea  surface  at  microwave  frequencies  can  be  calculated  by  assuming  that 
the  scattering  is  essentially  a  reflection  process  from  a  number  of 
bright  spots,  (stationary  phase  points)  characterized  bv  their  position 
and  height  above  the  mean  surface,  and  their  gaussian  curvature  g^. 
In  this  case,  the  scattering  cross-section  'q  of  a  portion  of  surface 
containing  N  bright  spots  is  given  by 


N 


I 

1 


1/2 


-jk(ri  +  r! ) 


(3-1) 


where  r^ ,  r]  are  the  distances  from  the  bright  spot  to  transmitter  and 
receiver  respectively  (see  Figure  3-1),  k  is  the  electromagnetic  wave 
number,  v  is  the  local  angle  of  incidence,  and  R  is  a  polarization 
dependent  reflection  coefficient.  Both  the  probability  distribution  for 
the  curvatures,  and  for  the  distribution  of  bright  spot  heights,  are 
known  for  the  sea  surface,  since  they  can  be  found  from  the  joint  height 
distribution  function  for  the  surface  [4].  Thus  it  is  possible,  using 
these  distributions,  to  produce  a  model  surface  which  is  composed  of 
a  number  of  randomly  distributed  spherical  or  ellipsoidal  sections 
whose  height  and  curvatures  are  appropriately  chosen  (see  Figure  3-1). 
This  model  should  provide  a  good  estimate  of  the  scattered  field  strength 
at  any  point.  Unfortunately,  the  problem  of  tracking  the  positions 
of  these  stationary  phase  points  as  both  transmitter  and  field  point 
move  is  a  formidable  one.  In  addition,  the  number  Na  of  points  per 
unit  area  of  surface  is  approximately  (for  back  scatter,  for  example) 


N  =  7-2  s2/(4Tth2)e'tan2v/s2 

d 


(3-2) 


for  a  plane  wave  incident  at  an  angle  v,  for  a  surface  with  rms  height 
h2,  and  rms  slope  s2.  A  sea  surface  area  of  even  a  few  thousand  square 
meters  would  possess  an  extremely  large  number  of  such  specular  points. 
For  these  reasons  (and  because  the  specular  point  model  does  not 
account  for  the  Bragg  scatter  component)  the  specular  point  model  was 
rejected  for  the  purposes  of  this  report. 
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TRANSMITTING  RECEIVING 


Figure  3-1.  The  specular  point  model  for  the  sea  surface. 


( i  i )  The  Explicit  Surface  Representation 

The  second  surface  model  to  be  considered  was  one  which  uses 
an  explicit  representation  of  the  surface  height  as  a  Fourier 
series.  For  a  square  area  of  side  L,  the  sea  surface  height  has  the 
expansion 


r.(xyt)  =  T  A  cos 


( 


where 

=  (2i/L)  (iZ+n2)1'  2 


is  the  mechanical  wave  number  and  wmn  is  the  angular  frequency  of  a 
component  traveling  wave  of  phase  ^mn-  Here  kmn  and  wmn  are  connected 
by  an  appropriate  dispersion  relation,  for  example  w^n  =  g  kmn  for 
gravity  waves.  However,  even  if  the  amplitudes  Amn  were  chosen  to  be 


5 


compatible  with  the  known  height  spectrum  of  the  ocean  (e.g.,  Amn 
B / ^nm  *or  J  Phillips  spectrum)  tins  representation  cannot  easily  he 
made  realistic.  This  is  because  the  coefficients  and  phase  angles 
must  be  chosen  in  such  a  way  as  to  reproduce  not  only  the  cruder 
statistical  properties  of  the  surface  (spectrum,  height  distribution) 
but  also  those  unique  features  of  the  real  sea  (cuspedness,  upwind 
skewed  slope  distribution,  etc.)  which  are  significant,  for  the 
scattering  problem.  Thus,  for  example,  the  problem  of  cuspedness 
might  be  approached  by  expanding  the  smaller  wave  number  components 
in  one  of  the  canonical  wave  profiles  ( i . e . ,  trochoids)  rather  than 
sinusoids.  However,  even  if  a  suitable  surface  representation  could 
be  found,  as  a  scattering  model  it  suffers  from  the  same  disadvantages 
as  the  specular  point  model.  This  is  that  the  direct  calculation  of 
the  scattering  from  such  a  surface  (e.g.,  by  integrating  "physical 
optics"  currents  over  the  surface)  would  require  impractical ly  long 
computation  times. 

(iii)  The  Flat  Patch  Model 

To  avoid  the  computational  difficulties  of  these  surface  models, 
it  was  decided  to  adopt  a  simpler  approach  to  the  surface  scattering 
problem.  The  model  chosen  is  the  "flat  patch"  model  in  which  the 
actual  sea  surface  is  replaced  by  a  flat  surface  divided  into  squares 
or  rectangles  (see  Figure  5-2).  The  scattering,  including  time  delay 
and  doppler  shift,  is  referred  to  a  conveniently  chosen  point  (usually 
the  center)  of  the  rectangle.  The  magnitude  of  the  scattered  field  in 
any  direction  is  then  determined  by  the  sice  of  the  rectangle  and  the 
square  root  of  the  bistatic  radar  cross-section  per  unit  area.  The 
ensemble  average  signature  produced  by  such  a  model  should  agree  fairly 
well  w’th  the  average  signatures  produced  by  the  specular  point  or 
surface  representation  models.  This  is  because  the  radar  cross- 
section  is  itself  derived  from  a  specular  point  model.  However,  to 
obtain  an  individual  signature  more  representative  of  an  actual 
signature,  it  may  be  desirable  to  associate  a  random  phase  shift  with 
the  return  from  each  patch. 

The  only  remaining  parameter  in  the  model  is  the  choice  of  patch 
sice.  In  practice,  the  area  of  the  patches  are  probably  most  conveniently 
chosen  empirically,  by  reducing  their  sice  and  increasing  their  number 
until  the  statistical  properties  of  the  return  show  no  significant 
change.  A  more  detailed  description  of  the  patch  model  is  given  in 
section  V. 

IV.  THE  "ONE-DIMENSIONAL"  SHIP-SEA  MODEL 

(i)  Introduction 

The  first  simplified  ship  sea  model  to  be  considered  was  an 
extension  of  the  Naval  Weapons  Center  Composite  Surface  (see  Figure  2-1), 
as  described  in  Reference  [1].  The  ship  is  represented  bv  a  flat  plate 
at  the  end  of  the  surface.  Even  this  simple  model  provides  examples  of 


the  four  principal  scattering  mechanisms  wh'ch  would  contribute  to  a 
radar  signature  in  a  real  encounter.  These  are  (i)  Bragq  type  bark¬ 
scattering  from  the  trapezoidal  structure,  which  represents  the 
"clutter"  signal  in  the  absence  of  the  ship  target;  (ii)  specular 
back  scatter  from  the  sinusoidal  structure,  which  contributes  to  the 
clutter  level  at  scattering  angles  near  normal  incidence;  (iii)  forward 
scatter  from  transmitter  to  surface  to  flat  plate  to  receiver,  which 
represents  the  double  bounce  corner  mechanism  (the  predominant 
interaction  between  sea  surface  and  target);  (iv)  the  field  diffracted 
by  the  edge  of  the  metal  plate,  which  represents  direct  diffraction  by 
the  target  alone. 

In  this  section  two  separate  approaches  to  calculating  the  sig¬ 
natures  of  this  one  dimensional  ship  sea  model  are  described.  In  the 
first,  a  physical  optics  technique  similar  to  that  of  References  1-3  is 
used.  In  the  second  approach,  the  methods  of  geometrical  optics  and 
the  geometrical  theory  of  diffraction  (GTD)  are  used  to  calculate  all 
mechanisms  except  the  Bragg  scatter  from  the  trapezoidal  structure.  The 
reason  for  this  dual  approach  is  that  the  first  method  (physical  optics) 
has  been  shown  in  Reference  1  to  provide  an  effective,  though  not 
particularly  rapid,  method  for  calculating  the  scattering  from  the 
composite  surface.  However,  in  calculating  the  signatures  for  more 
realistic  ship-sea  models,  it  is  expected  that  GTD  techniques  will  be 
used,  since  the  programs  are  faster,  and  since  GTD  is  known  to  be 
effective  in  computing  the  scattering  from  aircraft  and  similar  complex 
objects.  Thus  to  validate  the  GTD  approach  for  these  future  studies, 
it  was  considered  desirable  to  compare  the  GTD  results  with  physical 
optics  results  known  to  be  in  good  agreement  with  experiments. 

( i i )  The  Physical  Optics  Technique 

The  physical  optics  approximation  finds  the  fields  scattered  due  to 
an  incident  magnetic  field  H-j  by  using  the  approximate  surface  current 
distribution  Js  =  2r*i  x  Hi.  A  detailed  discussion  of  the  physical  optics 
fields  scattered  by  the  composite  surface  alone  is  given  in  Reference  [1]. 
The  application  of  physical  optics  to  the  ship-sea  model  of  Figure  2-1 
together  with  a  number  of  comparisons  with  experimental  signatures,  is 
given  in  detail  in  Reference  [5],  In  this  report  we  summarize  briefly 
the  principal  results  of  Reference  [5], 

Referring  to  Figure  4-1,  if  the  surface  with  generators  in  the 
x  direction  has  the  profile  z  =  h(y),  the  radar  trajectory  has  dive  angle 
antenna  depression  angle  ,■»,  and  the  range  to  the  mean  surface  along  the 
boresight  is  RQ,  then  the  receiver  voltage  (Reference  [5])  is 

2jV0  r  (Rocos(-.‘-z-)  +  ysinz-hcos,-) 

r  Jy  cos,*-[(h-Rocos»)^  +  (y+RQsi  nn 

|  (ycosa+hsir.iO^/Rp"  -  |(hcosr-ysinr)'7Rt^ 

*  e^  e  ^  dy  (4-1 ) 
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wi  th 


R1  =  (h'R0C0St')^  +  (y+R  sine)*' 

Rj2  =  (h-R^sinc)2  +  (y-RVos,;)2 

and  it  has  been  assumed  that  the  antenna  pattern  is  gauss i an  ’n  both 
planes  with  pattern  exp(-A  tan?n)  in  the  plane  defined  by  the  bore 
sight  and  the  x  axis,  and  pattern  exp(-B  tan?;.)  in  the  plane  defined 
by  the  bores ight  and  the  trajectory.  By  evaluating  the  integral  over 
the  surface  profile,  the  contribution  of  the  sinusoidal  (gravity  wave) 
structure  and  the  composite  surface  (sinusoidal  plus  trapezoidal 
structure)  can  be  found.  The  target  contribution  due  to  the  double 
bounce  corner  may  also  be  evaluated,  using  the  same  integral,  by 
exploiting  a  number  of  simplifications  introduced  by  image  theory. 

The  fourth  contribution  to  the  total  signature,  that  due  to  edge 
diffraction  must  be  calculated  by  GTl)  (see  section  IV-iii).  A  detailed 
derivation  of  the  form  of  these  four  components,  and  a  computer  code 
for  determining  the  resulting  signatures,  is  given  in  Reference  [5]. 

As  an  example  of  the  use  of  Equation  (4-1),  Figures  4-2  and  4-3 
show  the  four  components  of  the  total  signature  for  two  encounter 
geometries.  In  these  figures,  the  distance  scales  YY  represent  the 
y  co-ordinate  of  the  radar  with  respect  to  the  position  of  the  corner 
plate,  and  the  ordinate  represents  receiver  voltage  (in  dR)  with  respect 
to  a  nominal  reference  level.  In  each  figure  the  first  signature  is 
the  contribution  of  the  gravity  waves  alone,  and  the  second  signature 
is  the  backscatter  from  the  composite  surface  (sinusoid  plus  trapezoid 
surface).  Note  that  the  largo  lobe  towards  YY  -  0  is  an  artifact  of 
the  program  caused  by  terminating  the  integration  at  the  end  of  the 
surface.  The  third  and  fourth  signatures  in  Figures  4-2  and  4-3 
show  the  double  bounce  corner  and  the  edge  diffraction  terms  re¬ 
spectively.  The  fifth  signatures  show  the  predicted  total  signature 
(the  sum  of  the  four  components)  compared  with  signatures  measured 
at  the  Naval  Weapons  Center.  In  both  figures  the  absolute  level  of 
the  measured  signature  was  adjusted  so  that  the  measured  and  calculated 
levels  for  a  calibration  sphere  were  brought  into  correspondence.  Both 
signatures  include  the  effects  of  range  gate  and  receiver  processing 
(detector  and  low  pass  filter).  From  these  figures,  and  the  numerous 
signature  comparisons  given  in  Reference  [bl,  it  is  clear  that  the 
physical  optics  technique  can  give  good  agreement  with  experimental 
results  for  the  simplified  one  dimensional  ship-sea  model. 
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Figure  4-2.  Component  parts  of  the  voltage  signature  for  the  simplified 
ship  sea  model.  From  top  to  bottom  the  signatures  show  the 
gravity  wave  signature,  the  composite  surface  signature,  the 
double  bounce  corner,  the  edge  diffraction  contribution  and 
the  total  (composite  signature).  The  bottom  signature  is  com 
pared  with  an  experimental  run  for  the  same  geometry. 

(Note  sphere  return  at  YY  /v  -20)  Antenna  depression  angle 
,ii=60'’.  Dive  angleij'  =  0".  Plate  angle  n -0° . 


1  ^ 
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Figure  4-3.  Component  parts  of  the  voltage  signature  for  the  simplified 
ship  sea  model.  From  top  to  bottom  the  signatures  show  the 
gravity  wave  signature,  the  composite  surface  signature,  the 
double  bounce  corner,  the  edge  diffraction  contribution  and 
the  total  (composite  signature).  The  bottom  signature  is  com¬ 
pared  with  an  experimental  run  for  the  same  geometry. 

(Note  sphere  return  at  YY  <-»  -20)  Antenna  depression  angle 
.i*30".  Dive  angle  i/-=15".  Plate  angle  h=5". 
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( i  i  i )  introduction  to  the  GTO  Mode  l 

In  order  to  validate  the  computationally  efficient  methods  of 
geometrical  optics  and  GTO  for  application  to  more  complex  ship-sea 
models,  it  was  decided  to  apply  them  first  to  the  same  one-dimensional 
surface  and  target  (Figure  2-1)  discussed  in  the  previous  section. 

The  specific  geometry  analyzed  is  shown  in  Figure  4-4.  The  sea 
surface  is  represented  by  the  one  dimensional  sinusoidal  profile,  and 
the  simplified  ship  model  is  the  eight-foot  flat  plate  located  at  the 
end  of  the  sea  surface.  The  antenna  behaves  as  both  a  transmitter  and 
a  receiver  as  it  moves  along  the  indicated  antenna  path.  To  reduce 
the  complexity  of  the  computations,  only  the  case  of  the  u=90°  corner 
is  considered.  However,  it  is  not  difficult  to  modify  the  procedure  to 
handle  any  angle  u. 


Figure  4-4.  2-dimensional  geometry  of  the  simplified  ship-sea  model. 


The  major  scattering  mechanism  to  be  calculated  by  the  GTO  method 
is  the  double-reflected  field  between  the  sinusoidal  surface  and  the 
vertical  plate.  A  numerical  search  technique,  via  the  law  of  reflection, 
is  used  to  find  all  the  specular  points  on  the  sinusoidal  surface;  then 
all  the  rays  from  these  specular  points  to  the  flat  plate  and  back  to 
the  receiver  are  sunmed  to  obtain  the  scattered  field  at  the  receiver. 

In  addition,  the  diffracted  field  from  the  edge  of  the  plate  is  also 
calculated. 
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Since  the  vertical  plate  and  the  sinusoidal  surface  are 
cylindrical  structures,  and  both  the  antenna  boresight  and  the  antenna 
path  are  in  the  principal  plane  of  the  geometry  shown  in  Figure  4-4, 
a  two  dimensional  analysis  is  adequate.  In  addition,  since  the  antenna 
is  in  the  near  field  of  the  sinusoidal  surface  and  the  vertical  plate, 
a  near  field  analysis  is  required.  In  the  sections  that  follow,  a 
procedure  is  first  described  to  find  the  specular  points.  These  points 
are  then  used  to  calculate  the  double  bounce  (reflected)  fields,  and 
the  edge  diffracted  component. 

( i v )  Patterns  and  Incident  Fields  of  the  GTD  Model 

1.  Antenna  pattern  and  the  field  incident  on  the 
sinusoidal  surface 

The  antenna  pattern  is  assumed  to  have  a  15-degree  half-power 
beamwidth  with  a  gaussian  distributed  pattern  function.  Since  the  field 
points  are  in  the  near  field  of  the  antenna,  it  might  appear  that  a 
composite  source  (i.e.,  an  array  of  elements)  would  be  required;  how¬ 
ever  it  was  found  from  the  results  that  a  single  infinitesimally  small 
element  with  a  gaussian  pattern  function  was  sufficient.  This  simplifi¬ 
cation  saves  a  tremendous  amount  of  computer  time. 


The  gaussian  pattern  function  is  given  by 


f (0)  =  exp[-(ka  coso )  /4n ] 


(4-3) 


where 


k  =  2n/X 

a  =  constant  that  determines  the  beamwidth 
o  =  pattern  angle  defined  in  Figure  4-4. 


For  a  "horizontally"  polarized  antenna,  the  radiated  field  is  given 


E^e.s)  =  C,  •  f(o) 


(4-4) 


For  a  "vertically"  polarized  antenna  (the  E-field  vector  is  parallel 
to  the  plane  of  incidence),  the  radiated  field  is  given  by 
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E>,s) 


C,  si m'  f(i')  ■ 


i' 


(4-!0 


where  C\  and  Cv>  are  appropriate  constants,  and  s  is  distance  from  the 
source  to  the  field  point. 

2.  Law  of  reflection  and  the  numerical  search  technique 

In  order  to  find  the  field  reflected  off  the  sinusoidal  surface  in 
a  given  direction  d  as  shown  in  Figure  4-b.  the  reflection  point  must 
first  be  found.  The  point  of  reflection  must  satisfy  the  law  of  reflection 
which  states  that  the  angle  of  incidence  is  equal  to  the  angle  of  re¬ 
flection  that  is 

-  I  •  n  -  d  •  n  t4-t<) 

and  that  the  incident,  reflected  and  surface  normal  vectors  must  lie  in 
the  same  plane, 

1  x  n  -  d  x  n  ( 4  7 1 

I  -  unit  vector  from  source  to  reflection  point 

n  -  unit  surface  normal 

d  -  unit  vector  from  reflection  point  to  image  point. 


Figure  4-b.  Pouhle-rof led  ion  configuration. 
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These  equations  can  be  combined  as 

(n  x  I)(n  •  d)  +  (n  •  I)(n  x  d)  1  0 


(4-8) 


To  find  the  reflection  point,  the  following  expressions  are  needed 
(see  Figure  4-4): 


curved  surface  equation:  y(x)  =  h  cos 


f  (x_xo) 


(4-9) 


slope  of  surface:  ^  sin 


f  (x'xo) 


'  d  V  ' 

tangent  to  the  curve:  t  -  x  +  ^  y 


normal  to  the  curve:  n  -  -  ^  x  +  y 

(4-10a) 

-d.y/d\  x  +  v 

unit  normal:  n  =  - - 

(4- 10b) 

/ 1  +  (dy/dx)2 

incident  vector:  I  --  (x-x^)x  ♦  (y-.Vs )y 

(4-1  la ) 

reflection  vector:  d  =  (xf,-x)x  *  (y.-y).y 

(4-11 b) 

Substituting  ^guations  (4-10)  and  (4-11)  into  (4-8),  one  finds 


f  (nxIx+nyIy)(nxdx-nydx)t(nxdx+nydy)(nxIy-nvIx) 

=  t<vx>  dx  +  (y--vsH  [(*-*,.)  -  jj{-  (yr-y)]  ♦ 


+  ^x'xr)  dx  *  (yr-ynt(xs-x)  ♦  j-J  (ys-y)]  =  0 


(4-12) 


One  needs  to  find  all  the  reflection  points  on  the  sinusoidal 
surface  und?r  the  condition  that  the  point  (xr,y.-)  he  so  located  that 
the  vector  R  is  directed  toward  the  point  (xs,ys).  Since  the  wall  is 
a  flat  metal  surface,  the  vector  d  is  directed  straight  toward  the  image 
source  point  behind  the  wall  as  shown  in  Figure  4-5.  Hence  in  Equation 
(4-12)  one  should  replace  the  point  (xr,yr)  by  (-xs.ys). 

Since  Equation  (4-12)  is  very  complicated  and  very  time  consuming 
if  it  is  to  be  solved  repeatedly,  a  numerical  search  technique  (in¬ 
cremental  method)  is  utilized  to  find  all  the  specular  points.  The 
method  is  based  on  the  first  term  in  a  Taylor  series  expansion  of  the 
reflection-point  equations  (Equation  (4-12))  such  that  the  equations 
for  the  new  reflection  point  (fj+i)  is  given  in  terms  of  the  equation 
of  the  last  point  (fj)  by 


fj*i  *  fj  *  dfj 


where 


,>f  ,.f  jf( 

dfj  =  :>X  •'*$  *  ,.y  *  ;r«J  5,1 

S  S 


Because  fj+i  must  also  satisfy  the  reflection  equation,  it  implies 
that  fj  +  ]  -  11  or  fj  +  dfj  =  0  so  that 

>f,  .'f{ 


-  [  f  .  +  -  j-  x  +  ,<y  1 

J  Tx.  s  'y.  ■  sJ 


x  =  — 


'll 


(4-13) 


The  above  equation  says  that  as  the  source  moves  by  an  incremental 
distance  (0'xs.'''ys)  the  reflection  point  in  the  sinusoidal  surface  will 
move  by  an  incremental  distance  3x.  In  this  way  if  all  the  reflection 
points  are  once  found,  they  can  be  tracked  by  Equation  (4-13)  while 
the  source  changes  position  in  small  increments. 


Mn  the  event  that  xr  and  yr  would  re  replaced  by  the  coordinates  of 

the  image  of  the  source  as  imaged  by  the  reflecting  flat  plate. 
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( v)  The  Reflected  Field 


After  the  specular  points  have  been  found,  the  doubly  reflected 
field  for  each  specular  point  is  computed  using  the  geometry  illustrated 
in  Figure  4-6. 


IMAGE 


SOURCE  AND 


RECEIVER 


Figure  4-6.  Double-reflection  configuration. 


The  reflected  field  from  a  curved  surface  is  given  in  terms  of 
geometrical  optics  for  the  3-D  case  by  [6] 


Er(s)  =  E1 (Qp) 


r  r 


"1  °2 


-jks 


(r  i+s)(t->s+s) 

tl 
ii 

the  point  of  reflection  Qr  and  pT,  ^  are  the  principal  radii  of 


(4-1 


where  s  is  the  distance  traveled  by  the  reflected  ray.  In  this  case 
s  *  D2  +  D3  (see  Figure  4-6)  EMQr)  is  the  incident  electric  field  at 


1,2 


COSO 


r  .  2  .2 
sin  e,,  sin  o-j 

t 

[ 7 

r  .  2  .2-| 

sin  O2  sin  0^ 

R1  +  R2 

J 

cos2o1 

R1  +  R2  J 

where  s'  -  D1  is  the  distance  between  source  and  the  reflected  point 
o1  is  the  angle  of  incidence  with  coso1  =  -n  •  I  =  n  •  d 


t<,  is  the  angle  between  the  incident  ray  and  the  vector  tangent 
to  the  surface  in  the  x-y  plane  at  the  reflection  point 

is  the  angle  between  the  incident  ray  and  the  vector  tangent 
to  the  surface  in  the  z-plane 

and  R2  are  radii  of  curvature  of  the  curved  surface. 

For  the  two  dimensional  case,  R^  =-"  and  =  90°.  Thus  one  finds  that 

_  s'  R,  coso1 

0,  =  - j  and  0.  =  s' 

2s  '+R.|  coso1 


R  .  atr’ixLL1 

1  |y“(x)| 


1.5 


sign[y(x)] 


(4-15) 


R  is  the  dyadic  reflection  coefficient  such  that 


fi 


R  •< 


for  vertically  polarized  incident  wave  (electric 
vector  parallel  to  the  plane  of  incidence) 


1-1  for  horizontally  polarized  incident  wave  (electric 
v  vector  perpendicular  to  the  plane  of  incidence) 


It  is  necessary  to  sum  the  fields  calculated  by  Equation  (4-14)  from 
all  the  specular  points  for  any  instantaneous  source  location  to  form  the 
total  double-bounce  reflected  field.  All  fields  that  are  reflected  back 
to  the  receiver  are  multiplied  by  the  appropriate  antenna  pattern  function, 
giving 


received 


N  ,  /  of  •  D1 

l  E  (Qr)  1 - 

n=l  n  K  J(pj+D2+03)(D1+D2+D3) 


e-jk(D2+D3) 


f(0n) 


(4-16) 


The  field  traveling  in  the  reverse  direction  will  yield  the  same 
amount  of  energy  as  that  shown  in  Figure  4-6,  by  the  theorem  of  reciprocity, 
so  that  the  total  received  field  due  to  double  reflection  is  twice  that 
calculated  by  the  above  equation. 
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( vi )  Diffracted  Fields 

Any  ray  that  reaches  the  upper  tip  of  the  vertical  wall  will  excite 
a  set  of  edge  diffracted  rays  that  radiate  into  all  directions.  In  the 
present  analysis  there  are  two  major  types  of  diffracted  rays.  The  first 
type,  termed  a  reflected-di ffracted  ray,  is  associated  with  a  ray  reflected 
by  the  sinusoidal  surface  and  diffracted  by  the  tip  of  the  plate.  The 
second  type,  called  a  diffracted  ray,  is  associated  with  a  ray  that  comes 
directly  from  the  antenna  and  is  subsequently  diffracted.  The  6TD  is 
applied  here  to  analyze  these  two  terms.  In  order  to  define  a  number  of 
parameters  and  functions  used  in  the  GTD,  it  is  essential  to  discuss 
first  the  three-dimensional  wedge  diffraction  problem  (see  Figure  4-7). 

The  source,  which  radiates  the  field  E^s),  is  located  at  point  s'. 

It  can  be  an  arbitrary  electric  or  magnetic  source  causing  plane, 
cylindrical,  or  spherical  wave  incidence  on  the  edge.  The  diffracted 
vector  field  E^(s)  at  the  observation  point  s  can  be  written  in  terms  of 
a  dyadic  diffraction  coefficient  [7], 

Ed(s)  =  Ei(QE)  •  D  A(s)e'jks  (4-17) 

where 

D  =  D  D. 

o  o  s  h 

and  s  is  distance  between  diffraction  point  and  receiver.  The  quantity 
A(s)  is  a  ray  divergence  factor  given  by 


for  spherical  wave  incidence  where  s'  is  distance  between  source  and 
diffraction  point.  The  edge  diffracts  a  cone  of  rays  such  that  Scrfyj 
(see  Figure  4-7)  so  that  the  ray-fixed  coordinate  system  is  centered  at 
the  point  of  diffraction  Qe;.  The  relationship  between  the  orthogonal 
unit  vectors  associated  with  these  coordinates  (?o,^,0;  '  , I )  are 

given  by 

I  =  x  <fr'  D  =  i'Q  x  $. 

The  diffraction  coefficients  for  soft  and  hard  boundary  conditions  are 
given  by 
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Figure  4-7.  Geometry  for  three-dimensional  wedge  diffraction 


o  L 


+  cot 


^  ~2n  ~  ^F[kLa'  ' )]  ;  F[kLa+ (:+*)] 


+  cot 


(izhliL 

V  2n 


F[kLa  (i+t'U 


and  the  upper  sign  applies  to  the  soft  boundary  condition.  Hence 


F(x)  =  2j|J~xle‘^x  [  e 1  d'. 

1  IJ~I 

is  the  transition  function  and  involves  a  Fresnel  integral.  Also 


(4-18) 


a~(x)  =  2  cos 


2  2n-rN~-x 


in  which  N  are  integers  which  most  nearly  satisfy  the  equations: 
2tnN+  -  x  =  i 


2irnN  -  x  =  -ii 


The  distance  parameter  (L)  is  given  in  general  for  a  straight  wedge  by  [6] 


s ( ) p i p 2  sin  r> 

^g(n]+s)(p2+s) 


(4-19) 


where  pe  is  the  radius  of  curvature  of  the  incident  wavefront  at  Qe 
taken  in  the  plane  containing  the  incident  ray  and  the  unit 
vector  tangent  to  the  edge  at  Q^, 

p  1  * p 2  are  Pr^nciPa^  »~ad i i  of  curvature  of  the  incident  wavefront. 
For  the  two-dimensional  reflected-diffracted  case  considered  here  (see 
Figure  4-8a)  the  following  specific  values  and  equations  are  used: 
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whoro  R|  is  the  radius  of  curvature  of  the  sinusoidal  surface  and  is 
given  in  ( qua t ion  (4-l!>) . 

The  total  rofloc tod-diffracted  field  is  then  the  summation  of 
fields  coming  from  all  the  specular  points: 


ld(s)  = 
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n  I 
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(4 -20) 


By  the  theorem  of  reciprocity,  the  ray  that  travels  in  the  reverse 
direction  as  shown  in  figure  4-Hb  will  prodine  the  sumo  voltage  at  the 
receiver  terminal  as  that  calculated  for  figure  4-Ha.  Hence  the 
received  field  is  twice  the  value  given  by  Iquation  (4-20).  All  re¬ 
ceived  fields  should  also  be  multiplied  by  the  antenna  pattern  function. 

The  diffracted  field,  as  shown  in  figure  4-l>  can  also  be  calculated 
by  liquations  (4-17),  (4-18),  and  (4-19).  The  appropriate  parameters  for 
this  case  (figure  4-9)  are 
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so  that  the  distance  parameter  becomes  I  1)3/9. 

(vii)  Results 

Before  showing  the  computer  calculated  results  of  the  scattered 
field,  it  is  enlightening  to  examine  the  location  of  the  specular  point 
on  the  sinusoidal  surface  for  various  antenna  positions.  If  the  height 
of  the  surface  is  set  to  zero,  there  is  only  one  reflection  point  which 
is  located  at  the  intersection  of  the  flat,  plate  and  surface  (see 
figure  4-10).  This  solution  is  independent  of  the  source  position. 
Note,  however,  that  if  the  corner  angle  is  not  a  right  angle  (u/‘}0), 
the  reflection  point,  will  not  be  at  the  intersection,  and  will  depend 
on  the  source  position. 
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Figure  4-9. 


Direct  diffracted  field. 


Figure  4-10.  Reflection  point  for  u=90". 


As  the  amplitude  of  the  sinusoidal  surface  increases,  the  number  of 
reflection  points  also  increases.  To  illustrate  this  effect,  a  set  of 
specular  point  results  are  shown  in  figure  4-11  as  the  antenna  moves 
along  the  trajectory  shown  in  Figure  4-4.  It  is  observed  (Figure  4 - 1 1 c  > 
that  the  number  of  reflection  points  is  reduced  when  the  source  is  moved 
closer  to  the  vertical  plate. 

For  the  configuration  shown  in  Figure  4-4  with  ■<  60",  horizontal 
polarization,  and  a  lb-degree  antenna  beam,  the  received  signal  wave¬ 
forms  are  shown  in  Figure  4-12.  Figure  4-12a  shows  the  signature  of 
the  corner,  i.e.,  the  double  bounce  reflection  alone.  The  signature  of 
the  double  bounce  reflection  plus  all  diffraction  terms  is  shown  in 
Figure  4-12b.  Figure  4-12c  gives  the  measured  data  from  F.neounter 
Simulation  Laboratory  (LSI.).  The  difference  between  Figures  4-12b  and 
4-1 2c  in  the  region  0- NT- 60  is  due  primarily  to  back  scatter  from  the 
small  scale  structure  of  the  surface. 

The  GTD  received  signal  waveform  of  the  model  shown  in  Figure  4-4 
with  angle  a  equal  to  30"  is  shown  in  figure  4-13.  In  this  figure, 
the  calculated  result  is  again  compared  with  the  measured  waveform 
taken  at  the  Encounter  Simulation  Laboratory. 

r 

All  the  above  calculations artd  measurements  have  included  a  range- 
gate  function.  This  range  r*a Fo^ requires  that  each  contribution  to  the 
received  field  be  multiplied  by  a  function  of  the  distance  traveled  by 
the  field  between  the  source  and  the  receiver.  That  is,  each  term  in 
Equations  (4-16)  and  (4-20)  is  multiplied  by  the  range  gate  function. 

Two  types  of  fields  that  are  back-scattered  directly  from  the  sinusoidal 
surface  have  not  been  included  in  the  GTD  results.  One  is  the  field 
directly  reflected  by  the  sinusoidal  surface  alone;  however,  it  is  not 
difficult  to  include  this  term  if  it  is  desired  to  show  the  gravity 
wave  contribution  to  the  signature.  For  example.  Figure  4-14  shows 
the  signature  for  a  60  when  the  geometrical  optics  (bright  spot) 
contributions  from  the  large  scale  gravity  waves  are  included.  The 
other  is  the  back  scattering  due  to  the  small  scale  structure  (Bragg 
scatter)  on  the  surface.  The  GT0  method,  at  the  present  stage,  is  nob, 
able  to  analyze  the  Bragg-type  scattering.  ;  ..V 

The  GT0  model  is  employed  here  only  to  verify  its  validity  by 
comparison  with  the  physical  optics  model  which  was  discussed  in 
section  IV-ii.  By  comparing  a  number  of  signatures  calculated  by 
each  method,  it  is  clear  that  excellent  agreement  is  obtained  between 
the  two  approaches  for  the  interaction  field  waveforms  (i.e.,  the 
doubly  reflected,  reflect-diffracted,  and  diffracted). 
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Figure  4-13.  6TD  signature  for  a=30".  (a)  calculated  signature  with 

no  diffraction,  (b)  calculated  signature  with  diffraction, 
(c)  measured  data  from  ESL  (Run  No.  36).  Dive  angle  i]’=0. 
Corner  angle  n=  0 . 
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figure  4-14.  GTD  calculation  of  signature  including  the  direct 

reflection  from  the  sinusoidal  surface  for  the  same 
geometry  as  Figure  4-12. 


V.  THE  3-DIMENSIONAL  SIMPLIFIED  SHIP  SEA  MODEL 

( i )  Introduction 

In  this  section  we  begin  the  consideration  of  a  more  realistic  ship- 
sea  model .  Although  at  the  present  time  the  structure  used  to  represent 
the  ship  is  relatively  simple,  the  model  has  the  potential  of  handling 
considerably  more  complex  structures  without  any  fundamental  change  in 
the  computational  process. 

Here  the  radar  system  is  represented  by  a  transmi t/ receive  antenna 
of  specified  pattern  which  moves  along  a  specified  trajectory.  The  ship 
is  modelled  by  a  single  elliptic  cylinder  (see  Figure  5-1),  and  the  sea 
is  represented  by  the  flat  patch  model,  that  is,  by  a  flat  surface  with 
specified  bi-static  scattering  behavior.  Because  our  interest  is  in 
radar  wavelengths  much  smaller  than  the  characteristic  dimensions  of 
the  ship,  and  in  trajectories  which  transit  at  relatively  short  range, 
the  problem  is  formulated  here  in  terms  of  the  geometrical  theory  of 
diffraction.  The  basic  geometry  is  shown  in  Figure  5-2,  in  which  the 
antenna  is  traveling  in  an  arbitrary  direction  A  with  its  main  beam  axis 
(bore-sight)  pointed  in  the  direction  B  and  is  illuminating  the  sea 
surface  over  an  effective  area  of  W  x  L.  This  effective  area,  which 
is  dependent  on  the  transmitting  antenna  pattern,  is  broken  up  into 
smaller  sub-areas.  The  scattering  from  each  sub-area  is  entered  into 
the  GTD  program  as  a  polarization  dependent  scattering  matrix  with 
scattering  coefficients  referred  to  the  center  of  the  sub-area.  The 
scattered  field  from  each  sub-area  of  the  sea,  then,  appears  to  be 
generated  by  a  new  source  which  in  turn  illuminates  the  metallic  cylinder, 
and  from  which  it  is  reflected  to  the  receiver.  The  first  objective  is 
to  compute  this  doubly  reflected  field  (the  ship-sea  interaction  term), 
because  it  represents  the  major  contribution  to  the  overall  target 
signature.  Since  the  elliptical  cylinder  is  a  three  dimensional  structure, 
a  3-D  analysis  is  necessary.  In  addition,  a  near-field  analysis  is  still 
required  since  the  antenna  may  be  in  the  near  field  of  the  cylinder  at 
the  ranges  of  interest. 
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Figure  5-1.  Simplified  ship  model  at  sea. 


( i i )  Antenna  Pattern  and_  Field  Incident  Upon  the  Sea  Surface 

The  antenna  is  taken  to  be  a  plane  rectangular  sheet  of  current  with 
coordinate  system  as  shown  in  Figure  5-3.  Although  a  near-field  analysis 
is  required  to  specify  the  range  dependence  of  the  radiated  field,  it 
has  been  found  that  the  angular  dependence  can  be  approximated  simply 
by  the  far  field  pattern  function  of  the  source. 

The  radiated  far  zone  electric  fields  in  spherical  coordinates  are 
given  by  [8] 
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One  must  note  that  the  x-y-z  coordinate  system  specified  in  Figure  5-3  is 
independent  of  that  specified  in  section  IV. 
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Although  the  aperture  distribution  and  corresponding  antenna 
pattern  may  have  any  form,  for  convenience  we  have  assumed  that  the 
aperture  current  distribution  is  z-directed  and  gaussian,  given  by  [8] 


Jz(x',z') 


where  a  is  a  constant  parameter  which  determines  the  beam  width. 
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Since  the  gaussian  source  function  is  diminishingly  small  outside 
the  limits  of  integration,  the  integrals  can  be  approximated  by  letting 
the  limits  of  integration  go  to  infinity,  with  the  result  that 
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(5-2) 


Ef  =  0- 


For  an  x-directed,  gaussian  distributed  current,  the  electric  fields 
are  given  by 

-  (cos2o+sin20cos2f ) 


Eo  =  e  ^r  cosocos*  e  4,1 


( 5  -  3a ) 


•  ..  -  (cos2o+sin2ocos2$) 

E*  ’  IS  ^  s^  «  4" 


(5-3b) 


Since  the  antenna  generates  a  15-degree  broadside  (o=90°,  *=90°)  beam, 
and  cosO-cosf  is  vanishingly  small  near  the  broadside  direction,  E ^  is 
negligible  for  the  x  polarized  source  distribution. 

To  conclude,  the  fields  that  are  incident  upon  the  sea  surface  are 
given  by  Equations  (5-2)  and  (5-3b)  for  a  vertically  (Jz)  and  horizontally 
(Jx)  polarized  antenna,  respectively. 

( i i i )  Sea-surface  Scattering  Characteristic 

The  scattering  characteristics  of  each  sub-area,  utilizing  the 
flat  patch  model,  will  be  determined  by  the  bistatic  radar  cross-section 
of  the  ocean.  The  complete  bi -static  cross-section  of  the  ocean,  over 
all  possible  bi-static  angles,  has  never  been  measured,  so  that  for  the 
purposes  of  this  report  a  cross-section  based  on  theoretical  considerations 
must  be  utilized.  Fortunately,  because  of  the  geometrical  configuration 
of  the  problem,  the  significant  ray  paths  (from  antenna  to  sea  surface 
to  ship  to  antenna)  all  involve  predominantly  forward  scatter.  At 
microwave  frequencies,  the  specular  component  of  the  total  scattering 
cross-section  of  the  ocean  is  determined  by  a  reflection  (i.e.,  an 
optical  scattering)  process  controlled  by  the  slope  distribution  of  the 
surface.  Thus,  referring  to  the  geometry  of  Figure  5-4,  where  the 
incident  ray  lies  in  the  y’-z1  plane,  the  radar  cross-section  per  unit 
area  is  given  by 
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Figure  5-4.  Geometry  of  bistatic  cross-section. 
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where  Rij  is  a  polarization  dependent  Fresnel  coefficient,  and  p(tOdt' 
is  the  probability  that  the  surface  normal  n  makes  an  angle  between  i> 
and  i^+dt'  with  the  vertical.  This  equation  determines  the  bistatic 
cross-section  for  a  wave  incident  from  the  direction  defined  by  Oj , 
into  the  direction  defined  by  os  and  ^s.  These  two  directions  then 
define  the  direction  n  of  the  surface  normal  required  for  specular 
reflection  from  oi  into  es^s,  and  the  local  angle  of  incidence  v.  The 
reflection  process  also  transforms  the  polarization  of  the  incident 
wave.  However  it  is  more  convenient  to  specify  the  polarization 
transformation  in  tenns  of  a  specific  transmitter  state  with  unit 
polarization  vector  a  (which  may  be  complex v  a  •  a"*  =  1 )  and  a 
specific  received  polarization  state  b,  (B-b*=l).  In  that  case 
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R(v)  =  Rl(v)(a-s)(b-s)  +  R.,(a)(a-p.)(b-ps) 

where 


Ra(v)  =  (cosv-E)/(cosv+E) 

R„(v)  =  (>  cosv-E)/(t  cosv+E) 

E  -  (.  -sin2,)1/2 

.  =  complex  dielectric  constant  of  sea  water 
(£.jXk  )/sin  2v 
s  x  k^ 
s  x  ks 


The  angles  v  and  are  found  from 

cos  2v  =  coso.cos*'  -  sine  .sin-.1,  cos.*- 
is  i  s  s 

cos  e  =  (cos-'s+cose .  )/2  cosv 
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The  probability  density  p(e)  for  the  surface  slope  distribution  depends 
on  a  number  of  oceanographic  variables,  but  for  the  purposes  of  this 
report  is  assumed  to  have  the  simple  isotropic  form 

pU)  =  exp(-tan2;Vs2(l+2s2) )  (5-9 

ns 


s  =  mean  square  surface  slope. 

Thus,  apart  from  an  unknown  phase  the  effect  of  a  patch  sea  is  that  of 
a  polarization  and  amplitude  transforming  reflection  of  the  incident 
ray  into  a  ray  in  the  "reflected"  direction. 

( i  v )  Numerical  Search  Technique  for  Po i_nts 
of  Reflection 

Once  the  "reflecting"  properties  of  the  patch  of  ocean  are  known, 
we  may  return  to  the  problem  of  calculating  the  double-bounce  con¬ 
tribution  of  a  particular  patch.  From  any  one  of  the  sub-areas  to  the 
receiver  via  the  reflection  from  the  elliptical  cylinder,  there  must 
exist  one  reflection  point  on  the  cylinder  except  when  the  receiver 
is  in  the  shadow  region  of  the  cylinder.  To  find  this  reflection  point, 
the  laws  of  reflection,  as  defined  by  Equations  (4-6)  and  (4-7)  must 


be  satisfied.  When  these  laws  are  applied  in  a  sy •> lema t i c  manner,  it 
is  found  [9]  that  the  position  of  the  reflection  point  can  he  recovered 
from  a  single  unknown  parameter  a  (to  he  defined  later)  which  satisfies 
a  sixth  order  polynomial, 

+ct.<'  +c«  <^+c  h:.,./  +c,.<  * t.  o  (h  - 1 o) 

6  1)  4  .1  i  1  o 

The  coefficients  C0,C],  etc.  are  functions  of  the  source  location, 
receiver  location,  and  cylinder  parameters.  Since  the  derivation  of 
this  polynomial  is  quite  a  laborious  process,  and  it  has  been  derived  in 
the  reference  given  above,  it  will  not  be  derived  again  here. 


The  polynomial  technique  is  very  time  consuming  if  it  must  be  carried 
out  for  every  sub-area  at  each  point  of  the  trajectory.  Hence  it  is  used 
only  for  initial  acquisition  of  the  reflection  point  for  the  first  sub- 
area  of  the  Nth  row  (A^]  )  as  shown  in  figure  This  reflection  point 

is  then  tracked  by  a  very  fast  numerical  search  technique  through  many 
small  increments  (width=.\v )  to  the  new  reflection  point  for  sub-area 
This  numerical  search  technique  is  thereafter  applied  to  find  all  the 
reflection  points  for  the  sub-areas  in  the  same  row  (sub-areas  Aqj,  A^/j , 
etc.),  and  this  procedure  is  repeated  for  each  row. 

The  polynomial  method  can  be  considered  as  an  exact  solution  (except 
for  a  numerical  truncation  error).  In  the  numerical  search  technique,  a 
new  reflection  point  is  found  from  the  previously  determined  reflection 
point;  therefore  this  approach  can  have  an  accumulation  of  errors.  In 
order  to  control  this  accumulation  of  errors,  an  error  bound  is  defined 
such  that  if  it  is  exceeded  by  the  accumulated  error,  the  number  of 
increments  is  increased  (.u  being  reduced)  until  the  error  hound  is 
satisfied.  The  numerical  search  technique  is  briefly  discussed  below. 

We  first  define  (see  figure  5-2)  position  vectors  s,  o,  and  r,  where 
the  vector  s  is  the  position  of  the  center  of  sub-area,  vector  o  is  the 
position  of  observation  point  which  is  fixed  for  any  instantaneous 
antenna  location,  and  vector  r  is  the  position  of  reflection  point  on 
the  cylinder.  Thus, 

s  :  xs  x  +  ys.y 

o  =  xxtyy  +  zz 

0  O'  o 

r  V  f  yry  f  V’ 

A  parameter  v  is  introduced  to  define  the  reflection  point  on  the 
elliptical  cylinder  such  that 
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where 
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n  =  nxx  +  =  b  c°s  v  z  +  a  sin  vy  -  •  normal  vector  (5-12a) 

t  =  !xx  +  V  +  lzz  =  ^  "  s> 

=  (a  cos  v-x$)x  +  (b  sin  v-y$)y  +  zrz  •••  incident  vector  (5-1 2b) 

1$  =  Rxx  +  R^y  +  R^z  =  (5-r 

=  (x0'a  cos  v^x  +  <Vb  S1n  +  ^zo"zr^z  ■’*  reflection  vector 

(5-1 2c ) 

After  performing  the  dot  and  cross  products.  Equation  (5-11)  can  be 
separated  into  two  equations: 

f(v.xs,ys)  -  (\I/"yIy)(nxRy-nyRx)  ♦ 

*  (nxVVy><nxly-nylx)  "  0  (5-'3a> 


9(».Xs,ys.zr)  =  ("xV"y'y>"z  *  ‘"xVW'z  *  0 


(5-1 3b) 


By  substituting  Equations  ( 5-1 2b )  and  ( 5-1 2c )  into  Equation  (5-1 3b) .  one 
finds  that 
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(5-14) 


which  is  a  function  of  v,  x  and  y  .  For  Equation  ( 5-1 3a ) ,  one  can  apply 
a  technique  similar  to  that  discussed  in  section  4  to  obtain 


fJtl  =  fj  ♦  ^  *  0 


where 
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Hence, 
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(5-15) 


which  is  a  function  of  the  old  reflection  point  (v,xs,ys)  and  the  in¬ 
cremental  distances  (Axs,Ays).  The  new  reflection  point  in  the  x-y 
plane  is  now  defined  by  Vj+i  =  vj+Av  and  can  be  substituted  along  with 
the  new  source  location  (xs+Axs,Ays+  ys)  into  Equation  (5-14)  to  yield 
the  new  reflection  point  in  the  z-coordinate.  After  the  new  point  (j+1) 
has  been  determined,  it  can  be  substituted  into  Equations  (5-15)  and 
(5-14)  again  to  find  the  next  new  point  (j+2).  This  process  continues 
until  it  reaches  the  center  of  the  adjacent  sub-area  from  which  an 
additional  reflected  field  is  then  calculated. 


( v )  Reflection  from  the  Elliptical  Cylinder 


Once  the  field  scattered  from  the  center  of  each  sub-area  has  been 
calculated  and  the  corresponding  reflection  point  on  the  cylinder  located, 
the  reflected  field  incident  on  the  receiver  is  readily  computed  using 
the  geometry  illustrated  in  Figure  5-6.  The  scattered  field  for  each 
sub-area,  as  described  In  section  V-iii,  Is  treated  as  the  field  Incident 
on  the  cylinder.  The  reflected  field  is  given  by 
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where 

^(Qr)  is  the  'incident'  field  at  the  reflection  point  Qr, 

S  is  the  distance  between  reflection  point  and  observation  point,  and 


Oppp  are  principal  radii  of  curvature  of  the  reflected  wavefront 
c  at  the  point  of  reflection,  namely 
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where 


S'  is  the  distance  between  source  and  reflection  point, 
o1  is  the  angle  of  incidence  with  coso1  =  n-R  =  -n*i. 
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Figure  5-6.  Geometry  for  field  reflected  from  cylinder. 


o,  i§  the  angle  between  the  incident  ray  and  the  tangent  vector 
(e])  to  the  cylinder  in  the  x-y  plane  at  the  reflection  point, 
(see  Figure  5-7) 

is  the  angle  between  the  incident  ray  and  the  tangent  vector 
L  ( ez )  in  the  z-plane  (see  Figure  5-7) 

R]  and  R2  are  the  radii  of  curvature  of  the  elliptical  cylinder, 


;b2cos2v+a2sin2v)3/2 


Using  the  above  expressions,  one  finds  that 


2  sin202 

R^coso1 


R  is  the  dyadic  reflection  coefficient  given  by 


where  ei  is  unit  vector  perpendicular  to  the  plane  of  incidence  (see 
Figure  5-7)  and  e,i ,  e^  are  unit  vectors  parallel  to  the  plane  of  incidence 
and  reflection,  respectively. 


RECEIVER 


Figure  5-7.  Reflection  coordinate  transformation  system 
on  the  cylinder. 


The  incident  vector  I  in. Figure  5-7  can  be  transformed  to  the 
orthogonal  coordinate  system  n,  e] ,  e?  as 

'  i '  j  ^  ' 

I  =  -coso  n  -  si  no  cosuj  e^  -  si  no  sinw  e-|  (5-17) 

where  u>  is  the  angle  between  e2  and  the  projection  of  I  onto  the  plane 
containing  e]  and  e^  such  that 

1  •  e^  =  -si no1  cos«  . 


Since  it  is  also  true  that  I  •  =-cos02>  we  have 


COSu)  = 


The  unit  vector  e]  is  found  by  the  orthogonal  relationship: 


e,  =  e?  x  — —  =  z  x  — — 

1  2  |S|  |S| 

where  n  can  be  obtained  from  Equation  (5-1 2a) .  Hence 
el  =  elxx  +  el/ 


a  sin  v _ 

2  2  2  2 
b  cos  v+a  sin  v 


and  e 


b  cos  v _ 

ly  2  2  2  2 

J  b  cos  v+a  sin  v 


From  Figure  5-7  one  can  easily  transform  the  vectors  used  in  the  dyadic 
reflection  coefficient  into  the  x-y-z  coordinate  system, 


-  sinw  +  cosw  e^ 

-  sinw  z  +  cosui  •  elx  *  +  cosw  •  e1yy 


(5-1 8a) 


=  I  x  e. 


( 5-1 8b) 


e„  =  R  x  eA 


(5- 18c) 


(vi )  Summary  of  Solutions  and  Results 

For  each  sub-area  of  the  sea  surface  scattering  matrix,  the  field 
incident  upon  it  is  calculated  from  the  antenna's  radiation  functions 
as  expressed  in  Equations  (5-2)  and  (5-3b).  The  scattered  field  from 
this  sub-area  is  then  obtained  by  multiplying  the  above  incident  field 
by  a  scattering  coefficient  which  is  a  function  of  the  angle  of  incidence, 
angle  of  reflection,  sea  surface  characteristics,  polarization  of  the 
incident  field,  and  size  of  the  sub-area.  It  is  clear  that  the  anqle 
of  reflection  is  determined  from  the  location  of  the  reflection  point 
on  the  cylinder.  This  location  is  first  calculated  by  the  polynomial 
method  and  then  tracked  by  the  numerical  search  technique  as  the 
scattering  source  moves  from  one  sub-area  to  another.  Once  the 
scattered  field  from  the  sub-area  to  the  reflection  point  on  the  cylinder 
is  computed,  the  reflected  field  from  this  reflection  point  to  the 
receiver  is  readily  calculated  by  the  method  of  GTO  which  is  given  by 
Equation  (4-14).  The  final  solution,  i.e.,  the  total  received  field, 
is  then  a  summation  of  all  the  doubly  reflected  fields  for  the  sub- 
areas  that  are  illuminated  by  the  antenna. 

Presently,  a  computer  program  to  obtain  the  reflection  points  on 
the  elliptical  cylinder  has  been  completed  and  the  results  verified  in 
terms  of  hand  calculations.  As  an  example,  for  the  case  of  broadside 
incidence  on  the  cylinder,  as  shown  in  Figure  5-8,  the  locations  of  the 
reflection  points  on  the  cylinder  versus  the  sub-area  positions  are 
illustrated  in  Figures  5-9  and  5-10.  Figure  5-9  shows  the  reflection 
points  in  the  x-y  plane  as  measured  by  the  angle  v,  and  Figure  5-10 
shows  the  reflection  points  in  the  z-plane  as  measured  by  the  height  zr. 
For  the  case  of  off-broadside  incidence  (45°),  shown  in  Figure  5-11,  the 
reflection  points  versus  the  positions  of  the  sub-area  are  illustrated 
in  Figures  5-12  and  5-13.  The  curves  of  position  generated  for  broad¬ 
side  incidence  are  symmetric  with  respect  to  the  center  sub-area  as 
expected.  The  curves  in  Figures  5-9  and  5-10  are  nearly  straight  lines 
for  the  reason  that  the  broad  face  of  the  elliptical  cylinder  is  nearly 
flat.  The  abrupt  change  at  the  end  of  the  curves  in  Figures  5-12  and 
5-13  is  due  to  the  reason  that  the  reflection  points  have  reached  the 
narrow  end  of  the  elliptical  cylinder  where  the  surface  curvature 
varies  drastically. 

At  the  present  time  the  construction  of  a  computer  code  to  implement 
this  ship-sea  model  has  not  been  completed.  However,  a  number  of  pre¬ 
liminary  computations  have  been  made.  By  taking  into  account  only  the 
double  bounce  corner  contribution  (forward  scatter  from  transmitter  to 
ocean  to  ship  to  receiver),  a  number  of  typical  signatures  for  the 
geometry  illustrated  in  Figure  5-14  have  been  calculated.  In  Figures 
5-15  through  5-19,  typical  signatures  with  horizontal  trajectories  made 
with  a  15°  antenna  beam  at  a  depression  angle  of  60°  are  shown,  for  the 
case  of  the  elliptic  cylinder  target  with  semi-major  axis  of  25  length 
units  and  a  semi-minor  axis  of  5  length  units.  The  trajectory  altitude  is 
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40.5  length  units,  and  the  horizontal  range  (radar  to  target  center) 
varies  from  0  to  120  units.  The  parameter  in  the  signature  is  the 
approach  angle  which  varies  from  broadside  (SAH=Q")  to  head-on 
($AH=90°).  Although  the  general  form  of  the  signatures  is  quite 
similar,  being  strongly  influenced  by  the  assumed  antenna  pattern  and 
range  gate,  the  relative  received  voltage  varies  markedly  as  the 
approach  angle  changes. 

Although  the  full  implementation  of  the  model  remains  to  be  carried 
out,  particularly  the  introduction  of  a  back-scattering  cross-section 
from  each  patch  (to  model  the  no- target  clutter  level),  it  is  evident 
from  Figures  4-16  to  4-19  that  the  model  can  produce,  in  a  computationally 
efficient  way,  realistic  signatures  for  the  encounter  of  a  radar  system 
with  a  ship  at  sea. 
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Figure  5-9.  Reflection  points  (measured  in  angular  position  V)  versus 
the  sub-area  location  for  the  case  shown  in  Figure  5-8. 
The  sixth  row  in  the  scattering  matrix  is  used. 
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Figure  5-10.  Reflection  points  (measured  in  vertical  height  zr) 

versus  the  sub-area  location  for  the  same  case  shown 
by  Figure  5-9. 
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Figure  5-11.  Characteristic  view  of  elliptical  cylinder  and  position 
matrix  for  off-broadside  incidence  of  the  antenna.  The 
antenna  is  instantaneously  fixed  at  location  (x=39.0", 
y=39.0" ,  z=40. 5"0  and  a=3",  b=30". 
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Figure  5-14.  Configuration  of  elliptical  cylinder  at  sea. 

patch  size  =3x3.  number  of  patches  on  each  side 
of  the  Illuminated  area  =  17.  Vertical  polarization 
is  employed. 
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Figure  5-15.  Voltage  signature  for  the  model  of  Figure  5-14.  Horizontal 
trajectory  (AT=0)  at  altitude  40.5  units  (ZF=40.5).  Ship 
dimension  A=5,  B=25  units.  Antenna  depression  angle  ATA= 
60°.  Aspect  angle  SAH=0°. 


VI.  SUMMARY  AND  CONCLUSIONS 

In  this  report,  two  models  for  calculating  the  interaction  between 
a  radar  system  and  a  ship  at  sea,  both  adapted  to  short  range  encounters, 
have  been  presented.  The  first,  developed  to  parallel  a  series  of  ex¬ 
periments  at  the  Encounter  Simulation  Laboratory  of  the  Naval  Weapons 
Center,  utilized  a  one-dimensional  composite  surface  to  model  the  sea, 
and  a  flat  plate  to  model  the  ship.  Four  major  scattering  mechanisms 
were  identified,  and  signatures  calculated  using  first  a  physical  optics 
formulation,  and  second  a  computationally  more  efficient  formulation 
based  on  the  geometrical  theory  of  diffraction. 

Signatures  based  on  the  physical  optics  method  were  in  excellent 
agreement  with  measured  signatures,  and  in  both  calculated  and  measured 
results  the  four  components  of  the  total  signature  could  be  clearly 
identified.  The  GTD  signatures  did  not  include  the  Bragg  scatter 
component,  but  when  compared  with  the  physical  optics  signatures  with 
Bragg  scatter  omitted,  again  excellent  agreement  was  obtained.  Thus 
it  was  concluded  that  the  GTD  method  could  produce  accurate  results 
for  all  but  the  Bragg  back-scatter  mechanism. 

The  second  ship-sea  model  utilized  a  flat  patch  model  for  the 
ocean  surface,  with  the  scattering  from  each  patch  determined  by 
the  bistatic  radar  cross-section  of  the  ocean.  The  ship  was  modelled 
by  an  elliptic  cylinder.  The  major  scattering  mechanisms  for  this 
model  (predominantly  the  double  bounce  corner  mechanism)  have  been 
identified,  and  efficient  methods  for  calculating  the  scattered  fields 
(based  on  tracking  stationary  phase  points  on  the  ship  structure)  have 
been  developed.  Some  preliminary  signature  calculations,  showing  the 
effect  of  the  aspect  angle  of  the  ship  target,  have  beer  made. 

However,  to  produce  realistic  signatures  for  this  model,  two 
additional  features  must  be  programmed.  The  first  is  to  include 
edge  diffraction  from  the  edges  of  the  cylinder,  and  reflection  from 
the  top  surface  of  the  cylinder,  as  the  trajectory  passes  over  the 
target.  The  second  is  to  include  the  back-scattering  from  each  flat 
patch,  in  order  to  establish  the  clutter  level,  and  typical  clutter 
signatures,  in  the  part  of  the  trajectory  preceding  the  encounter. 

Once  these  components  are  included,  it  should  be  possible  to  obtain 
realistic  signatures  for  a  variety  of  encounter  geometries,  for  systems 
with  any  prescribed  antenna  patterns  and  polarizations,  and  receiver 
processing  technique. 

For  future  development  of  the  model  one  must  consider  the  use  of 
a  more  realistic  ship  model,  for  example  a  hull  formed  from  a  section 
of  an  elliptic  cone,  perhaps  with  a  roll  angle,  and  surmounted  by  a 
rectangular  or  elliptic  superstructure.  In  addition,  one  should  con¬ 
sider  a  more  realistic  sea  surface  (including,  e.g.,  a  systematic  bow 
wave  structure  as  well  as  the  random  sea  surface).  The  techniques 
developed  in  this  report  should  permit  the  calculation  of  signatures 
which  include  both  refinements  in  the  ship  model  and  the  sea  surface 
model  suggested  above. 
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